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Abstract 



We discuss a simple method to evaluate the QCD corrections to Ap. It assumes that 
the perturbative expansion in terms of MS parameters is meaningful and, unlike other 
studies, exploits significant available information concerning 0(a 2 ) corrections. This 
approach leads to an enhancement of ~ 26% relative to the conventional evaluation. 
fNJ ■ QCD corrections to the Z° — > bb amplitude are also considered. Implications for 

CO . electroweak physics are briefly discussed. 

On 



The 0(aa s ) corrections to the electroweak amplitudes Ap and Ar have been the 



o ,| object of detailed studies in the past [1]. These analyses have been recently extended to 
the contributions proportional to m\ in the Z° — > bb amplitude [2]. In a recent paper, 
Smith and Voloshin (S-V) have re-examined the 0(aa s ) corrections to Ap by a detailed 
investigation of the relevant Feynman diagrams involving the t-b isodoublet [3]. Because of 
the large mass of the top quark and the fact that Ap is evaluated at a much smaller scale, 
namely q 2 = 0, the conventional wisdom is that the Feynman integrals are dominated by 
euclidean gluon momenta k ~ rat [4]. In fact, in most of the published calculations, a s is 
evaluated at mt, with the cases a s {2mt) and a s (mt/2) sometimes considered to estimate 
the theoretical error. However, in the case of Ap, S-V find a significant sensitivity to 
scales ~ 0.15mt. Specifically, employing the V-scheme running coupling [5], carrying out 
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the integration over the euclidean gluon momentum and absorbing the result in a rescaling 
of a s , these authors find: 



(Ap) f =x t 



2d s (0.154m t ) 2 . 

— (7T + 3) 



9tt 

where the subscript / means fermionic contribution, 



(1) 



is the one-loop correction [6] and m t the on-shell mass. Eq. (1) is the usual perturbative 
result [1], except that the MS coupling a s is evaluated at 0.154m i7 rather than the scale 
m t conventionally employed in these calculations. The authors of Ref. [3] point out that 
a complete three-loop 0(aa 2 s ) calculation is needed to completely quantify this scale. 
However, as things presently stand, the shift in the a s scale from m t to 0.154mt is quite 
large, and implies a significant modification in the 0(aa s ) corrections. For instance, for 
a s (mz) = 0.118 and m t = 200 GeV (roughly the upper bound of the current indirect 
determination of mt), d s (0.154mt) /a s (mt) = 1.34, where we have employed a three-loop 
evaluation of a s . Thus, the S-V rescaling of a s leads to a 34% enhancement of the usual 
QCD contribution to (Ap)f. Furthermore, as the corrections of 0((a s /n) 2 ) involved in 
this rescaling have large coefficients ~ —20.5, one wonders whether the neglected 0(6r s ) 
corrections may also be sizable. 

In this communication we discuss a simple method to evaluate the QCD contributions 
to (Ap) f which, in contrast with other studies, exploits significant available information 
concerning 0{a s ) corrections. As we will see, it also assumes that the perturbative expan- 
sion of (Ap) f in terms of MS parameters is meaningful. We start with the observation, 
made by S-V, that the sensitivity to relatively low mass scales exhibited in Eq. (1) arises 
from the use of the on-shell mass m t in Eqs. (1) and (2). In fact, these authors traced such 
dependence to the contribution to (Ap)f of the mass counterterm associated with m t . In 
order to circumvent this problem, we propose to first express ( Ap)f in terms of m t (m t ), 
the running MS mass evaluated at the m t scale. One readily obtains 



(Ap)/ = 



i-|^-.) + c(|) a + 



(3) 



where m t is henceforth an abbreviation for rh t (m t ), a s is also evaluated at m t , C is the 
coefficient of the 0((a s /Tv) 2 ) contributions which have not been calculated so far, and 
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the ellipses represent higher order terms. As the answer is expressed in terms of m t , the 
corrections in Eq. (3) should not be sensitive to the relatively low mass scales found by 
S-V and, in fact, the dominant contributions are expected to involve euclidean momenta 
k > m t . Furthermore, we note that the coefficient of the 0(a s /iv) term is small, namely 
—0.1932. We now recall that the relation between m t and m t is accurately known [7]: 



T 1 =1 + — ^ + 



16.11 



1.04V 1-^ 

feV rut 



— I + 

TV 



(4) 



mt 3n 

where a s is again evaluated at m t . For 130 GeV < m t < 200 GeV, we approximate the 
small mi-dependent term within the square brackets by 0.04 and rewrite Eq. (4) as 



^ = 1 + ^ + [16.15 - 1.04^] ^ 

mt 3 tv \ tv ' 



(5) 



where n/ = 5 is the number of light quark flavors. Numerically, the coefficient of the 
0((6l s /tv) 2 ) term is quite large, namely 10.95, which may partly reflect the sensitivity to 
the relatively small mass scales discussed before. In order to improve the convergence of 
the expansion, we apply the BLM method [5]. Writing n/ = 33/2 — 3/3o/2, where (3q is the 
coefficient of the one- loop (3 function for ct s , one absorbs the (3q term in a shift in the a s 
scale, while —1.04 x 33/2 is combined with 16.15. Eq. (5) becomes then 



m t 3 tv 



1 - 0.758 



TV 



+ 



(6) 



where Q* = 0.0963mt. This procedure absorbs vacuum polarization contributions in the 
coupling constants. We see that Q* <C m t , a frequent feature of the BLM method, but at 
the same time the convergence of the series is much improved. For example, for m t = 200 
GeV and a s (m z ) = 0.118, we have a s (m t ) = 0.1055 and a s (Q*) = 0.1546, so that Eq. (5) 
gives 1+0.0448+0.0123 + - • • = 1.0571 + - • •, while Eq. (6) leads to 1+0.0656-0.0024 + - • • = 

1.0632 -| . We also note that the two-loop result from Eq. (6) is larger by 6 x 10~ 3 . As 

Eq. (6) contains a very small next to leading term, it gives essentially the same result as the 
method of fastest apparent convergence (FAC) [8]. The latter would absorb the complete 
0{a 2 s ) contribution in Eq. (5) into a rescaling of the 0(a s ) term. In the following, we 
employ Eq. (6). 

Combining Eqs. (3) and (5), we have 



(Ap) f = x t 



1-f ^ 2 ~9) + ^) : 
(m t /m t ) 2 



(7) 
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where a s in the numerator is evaluated at m t and m t /rh t is given by Eq. (6). For a very 
precise evaluation of Eq. (7), the determination of the constant C would be needed and this 
is not available at present. However, if we make the plausible assumption that the 0(a 2 ) 
term in Eq. (3) and the numerator of Eq. (7) is not larger than the corresponding 0{a 8 ) 
contribution, it is clear that the bulk of the QCD correction resides in the denominator of 
Eq. (7). In fact, this factor contains the dependence on the relatively small mass scales 
discussed above and, moreover, it clearly absorbs most of the 0(a 8 ) contribution. It is 
also interesting to note that if one expands Eqs. (1) and (7) in powers of a s (mt) up to 
terms of 0(a 2 ), one can match the second-order coefficients by choosing C = —4.46. This 
means that the denominator of Eq. (7) absorbs most of the 0(a 2 ) contribution of Eq. (1) 
(— 16.1(d s /7r) 2 out of — 20.5(a s /7r) 2 to be precise). The important theoretical point is that 
this large contribution, (m t /rh t ) 2 , has been evaluated with considerable accuracy (Eqs. (5) 
and (6)). On the basis of these arguments, we propose to calculate (Ap) / using Eq. (7) with 
a s evaluated at m t in the numerator, and m t /m t given by Eq. (6). As the constant C has 
not been determined, we treat the last term in the numerator of Eq. (7) as a measure of the 
theoretical error associated with unaccounted higher order corrections. Assuming that the 
perturbative expansion in Eq. (3) is meaningful, we fix \C\ so that the 0(a 2 ) term in that 
equation is not larger in absolute value than the 0(a s ) contribution. This leads to C ±6. 
With this choice, the coefficient of (a s /n) 2 in Eq. (3) and the numerator of Eq. (7) may be 
as large as 6/0.1932 ps 31 times the cofactor of a s /n. By way of comparison, in the case of 
T(T — > hadrons)/r(T — > p+p - ), which is often cited as an illustration of large higher order 
coefficients, the cofactor of the next to leading term in an expansion in powers of a s (Mx) /V 
is 9.1 relative to the leading contribution; in the case of r(rj c — > hadrons)/r(?7 c — > 77), 
an expansion in powers of a s {M r]c ) /n shows a 14.5 relative coefficient [5]. We may also 
attempt to estimate the vacuum polarization contribution to C by combining Eqs. (1,5,7): 
one finds C vac .poi. ~ 0.44,30. As in this case the relevant scale is k > m t , it is natural 
to consider six active quark flavors so that C vac . po i. ~ 3.1. (These "vacuum polarization" 
effects include gluonic contributions associated with the running of a s ). Another way of 
estimating C vac . po i. is to apply the BLM method and absorb this contribution in a rescaling 
of the 0(a s ) term in Eq. (3). Evaluating accurately the rescaled term one finds that, for 
m t = 200 GeV, the result is equivalent to using an effective value C va c.poi. ~ 2.2 in Eq. (3). 
We note that both estimates are positive, which is consistent with the notion that the 
dominant contributions to Eq. (3) involve k > m t . Yet another approach to estimate C is 
to apply the principle of minimum sensitivity (PMS) [9]. One expresses Eq. (3) in terms of 
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rht(fi) so that (Ap)f is proportional to mf (p){l + (a s (p)/Tv)[4£n(p/m t ) - (2/9)(n 2 - 9)]}. 
Evaluating this quantity at the stationary point p* , which according to the PMS is the 
optimal value, we find that, for m t = 200 GeV, the result is equivalent to C 2.7 in 
Eq. (3). Thus, we see that these estimates of C are roughly consistent and all comfortably 
lie within our error range. On the basis of the various arguments given above, it appears 
that our choice C = ±6 in the error analysis is quite conservative. Another possible source 
of theoretical uncertainty arises from the 0({a s (Q*)/n)' i ) contributions to Eq. (6). If, 
following the pattern of Eq. (6), the corresponding coefficient is ~ 1, such terms would 
give very small effects ~ ±2 x 10 -4 . If the coefficient is ~ 10, their effect would be 
~ ±2 x 10 -3 , three times smaller than the estimated error we have included. Using then 
a s = a s (m t ) = 0.1055, a s (Q*) = 0.1546 (values corresponding to m t = 200 GeV and 
«s(^iz) = 0.118), and C = ±6, Eq. (7) gives, for the QCD correction to {Ap)f. 



Eq. (8) is to be compared with the conventional result (Ap) f — x t = — 0.0960x t (Eq. (1) 
with a s (m t )) and the S-V proposal ( Ap) f — x t = — 0.129x t (Eq. (1) with d s (0.154m t )). 
Thus, Eq. (8) leads to a 26% enhancement of the QCD correction over the conventional 
answer, and a 5% error estimate, relative to the calculated QCD effect, due to the higher 
order corrections. It is somewhat smaller than the S-V result, but roughly consistent with 
it if one takes into account the estimated uncertainty. A more complete discussion of 
the error must include the effect introduced in Eq. (8) by the experimental uncertainty in 
a s {mz)- For 5a s (mz) = ±0.006, we find an additional error ±0.007 in Eq. (8). Combining 
the two errors in quadrature, the overall uncertainty in Eq. (8) becomes ±0.009, or about 
8% (if added linearly the total error would be ±0.013 or 11%). 

For mt = 200 GeV, the difference between Eq. (8) and the conventional evaluation 
leads to an additional contribution —3.1 x 10 -4 to ( Ap)/ and +1.1 x 10 -3 to Ar [10]. 
The latter corresponds to a decrease ~ 19 MeV in the predicted value of for given 
m t , or an increase Am t = ±2.8 GeV for given mw- On the other hand, the S-V proposal 
(Eq. (1)) implies somewhat larger shifts, namely Amw = —25 MeV and Am t = 3.7 GeV. 

An analogous argument can be applied to the term proportional to mf in the Z° —> bb 
amplitude. Including the 0(aa s ) correction, it is of the form [2] 



(Ap) f -x t 



x t [0.121± 0.006]. 



(8) 




m\ [1-(tt 2 -3) d s /37r + •••], 



(9) 
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where the superscript v means "vertex" and the subscript b refers to Z° — > bb. In terms 
of m 2 , this becomes rhf [1 + (11 — 7r 2 )d s /37r + D(a s /n) 2 + •••]. Thus we see that, once 
more, the coefficient of the 0(a s ) term is significantly reduced. The latter expression can 
be rewritten as 

2 K (v) 2 [l + (ll-7r 2 )fe+D(« s /7r) 2 ] 
m t Kl'=m t f—^ , 10 

where again a s in the numerator is evaluated at nit and nit I nit in the denominator is given 
by Eq. (6). For nit = 200 GeV, we choose D = ±11 (which corresponds to the situation in 
which the 0(a 2 ) and 0(a s ) terms in the numerator are equal in magnitude), and obtain, 
for the QCD correction, m 2 (K { b v) - 1) = -m 2 [0.104 ± 0.011]. This is to be compared with 
the conventional result — m 2 0.0769. So, in this case we find a 35% enhancement over the 
conventional answer and an 11% error estimate, relative to the calculated QCD effect, due 
to the higher order corrections. 

It should be pointed out that the use of the on-shell mass in Eq. (2), which is the 
usual choice, implicitly assumes that this is the relevant parameter to discuss top physics: 
its production, decay and the measurement of its mass [11]. For instance, if the running 
mass nit were found to be an appropriate parameter to analyze certain aspects of top 
physics, one could directly use Eq. (3). The idea of expressing (Ap)j- in terms of m 2 , with 
or without additional 0(a s ) corrections, and to employ Eq. (5) up to 0(a s /n) or higher 
to relate the two masses, has occurred to a number of authors [11,12]. Our strategy has 
been to separate the QCD corrections to (Ap)f into two parts: 1) The factor {mt/rht) 2 , 
where the corrections are large but known through 0(a 2 ). We have further shown how the 
convergence of this series can be significantly improved by applying the BLM method, and 
noted that essentially the same result is obtained in the FAC approach. Thus, the present 
framework employs significant information concerning 0(a 2 ) contributions. 2) The second 
part involves the corrections to (Ap)/, when expressed in terms of m\. Here the 0(a s ) 
contribution is small, while the 0(a 2 ) term has not been computed so far. By assuming 
that the latter is not larger than the former, so that perturbation theory may be viable, the 
magnitude of the 0(a 2 ) term in that expansion is bounded. We have also noted that this 
procedure allows for the possibility that the coefficient of the (d s /7r) 2 term may be as large 
as 31 times the cofactor of the a s /it contribution. In spite of this latitude, we obtain a 
rather small estimate for the theoretical error of the overall perturbative QCD correction. 
Analogous comments apply to the analysis of the QCD corrections to the m 2 terms in the 



6 



Z° — > bb amplitude. Our argument could go wrong in two ways: i) If the perturbative 
expansion in Eq. (3) and the analogous one for the Z° — > bb amplitude fail. This would 
be the case if the coefficients C and D are significantly larger than the values \C\ = 6 and 
\D\ = 11 we employed in the argument, or if, for special reasons, terms of nominal 0(a 3 ) 
or higher are sufficiently large [13]. The elucidation of these issues lies outside the scope 
of this paper. Here we have adopted the point of view that the perturbative expansion in 
Eq. (3) and the corresponding one for the Z° — > bb amplitude are meaningful, which seems 
very reasonable in view of the large mass scales we are considering. This point of view is 
supported by well-known theoretical arguments invoked in the derivation of QCD sum rules 
for heavy quarks [4,14]. If, in the future, additional significant contributions to Eq. (3) can 
be unambiguously shown to exist, they would have to be added to the perturbative effects 
we have discussed, ii) If the coefficient of the 0((a s (Q*)) 3 ) term in Eq. (6) is considerably 
larger than 10. In this regard, we point out that the nice convergence pattern of Eq. (6) is 
not suggestive of an explosive behavior in the 0((a s (Q*)) 3 ) contributions. Nonetheless, it 
is clear that these studies will be greatly strengthened when complete 0(aa 2 s ) calculations 
become available and the constants C and D are determined. They will also greatly benefit 
from the analysis of 0((a s (Q*)) 3 ) corrections to Eq. (6). Meanwhile, Eqs. (6), (7), and (10) 
provide, in our view, a suitable perturbative framework to evaluate the QCD corrections 
to (Ap) f and Kl v) . 
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